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Abstract 

We quantize the Reissner-Nordstrom black hole using an adaptation of Kuchaf 's 
canonical decomposition of the Kruskal extension of the Schwarzschild black hole. 
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The Wheeler-DeWitt equation turns into a functional Schroedinger equation in 
Gaussian time by coupling the gravitational field to a reference fluid or dust. The 
physical phase space of the theory is spanned by the mass, M, the charge, Q, the 

a ■ 

physical radius, R, the dust proper time, r, and their canonical momenta. The 
exact solutions of the functional Schroedinger equation imply that the difference in 
the areas of the outer and inner horizons is quantized in integer units. This agrees 
in spirit, but not precisely, with Bekenstein's proposal on the discrete horizon 
area spectrum of black holes. We also compute the entropy in the microcanoni- 
cal ensemble and show that the entropy of the Reissner-Nordstrom black hole is 
proportional to this quantized difference in horizon areas. 

PACS: 04.60.Ds, 04.70.Dy 



f Email: cvaz@ualg.pt. 

| Email: witten@physics.uc.edu 



1 



I. Introduction. 

Although the temperature of a black hole is exactly zero degrees Kelvin in 
classical general relativity, Bekenstein, in his 1972 thesis, 1 proposed that black holes 
have a temperature and entropy, and should be treated as thermodynamic systems. 
The temperature and entropy of the black hole are known from semi-classical 
arguments 2 ' 3 ' 4 to be fundamentally quantum mechanical in nature. Therefore, 
understanding their origins from a bona fide microcanonical ensemble of quantum 
states has come to be recognized as a challenge of considerable importance for an 
eventual theory of quantum gravity. 

The earliest attempt at a microscopic theory of black holes was also due to 
Bekenstein. The argument roughly goes as follows. 5 ' 6 Using the Christodoulou- 
Ruffini 7 process, it becomes clear that the horizon area operator of the black hole 
(in the case of multiple horizons, the area of the outer horizon) must be treated as 
an adiabatic invariant. Then, invoking the semi-classical Bohr-Sommerfeld quanti- 
zation rules, Bekenstein concluded that the horizon area operator admits a discrete, 
equally spaced spectrum, A n ~ n/ 2 , where l p is the Planck length, and proceeded 
to use this spectrum as the rationale for dividing the horizon into cells of unit area 
which get added one by one and which have the same, small number of states, say 
k. The result is an estimate of the density of microstates, f2 ps /c n , or the entropy 
of the black hole, S — lnfi w nlnk In k(A/lp). 

This paper is a development of earlier work, 8 ' 9 where the quantum states and 
the total entropy of the Schwarzschild black hole were recovered by combining the 
canonical reduction of spherical geometries by Kuchaf 10 with the coupling to an 
external reference fluid (or dust), originally proposed by Kuchar and Torre 11 and 
Kuchar and Brown. 12 The functional Schroedinger equation (in the dust proper 
time) that was obtained by this procedure described the more general problem 
of inhomogeneous dust collapse. It was simplified by holding the mass of the 
hole constant, independent of the spatial coordinate, and could then be easily 
solved throughout the Kruskal manifold. This gave precisely Bekenstein's area 
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quantization law. The coupling to dust may be thought of either as a way to 
impose coordinate conditions (ref. [11]), which is the point of view we take here, 
or as a realistic material medium (ref. [12]). 

Our goal in this paper is to extend this analysis to the charged, Reissner- 
Nordstrom black hole. We will see that, as in the Schwarzschild case, the quan- 
tization leads to a derivation of the statistical properties of the black hole. In 
particular, the entropy will turn out to be proportional to the difference between 
the outer and inner horizon areas and will be quantized in integer units. Thus, 
although we do not recover Bekenstein's area quantization, our result is in keeping 
with its spirit in as much as it is commensurate with an "area quantization" law. 
As the charge, Q, approaches zero the results will approach those obtained for 
the uncharged Schwarzschild black hole and in the limiting the black hole 

becomes extremal, the entropy will approach zero. 

The Reissner-Nordstrom solution is given in curvature coordinates by the line 
element 

ds 2 = F(R)dT 2 - F-\R)dR 2 - R 2 dtt 2 , (1.1) 

where dfl is the ordinary unit sphere, R is 
has the form 

F(R) = 1 - 
and the electromagnetic potential is 

A = ^dT. (1.3) 

The vector field d/dT is a Killing vector field of the metric. It is time-like in the 
interior (regions IV and V, see the Penrose diagram in figure 1) and in the exterior 
(regions I and II), but space-like in region III. 

An important feature of the maximal extension of this geometry is that the 
inner horizon is a Cauchy horizon for spatial sections, E (see figure 1). 



the physical radius, the coefficient F(R) 
2M Q 2 , x 

IT + W < L2) 
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Fig. 1: The extended Reissner- Nordstrom geometry 

What this means is that if data were given on an initial hypersurface, S, the 
Cauchy development will be able to predict only what occurs in regions I, II and 
III and not beyond the inner horizons, in regions IV, V of the diagram. Any event 
in regions IV and V of the space-time is influenced not simply by the given data 
and evolution but by additional data on the singularities themselves, which are 
impossible to control. Such a situation does not arise for the black hole (see figure 
2) where spatial sections are able to cover all of space-time until the singularity is 
reached and the Cauchy development is able to predict what occurs everywhere, 
once data is given on an initial hypersurface. 



r=0 




r=0 

Fig. 2: The extended Schwarzschild geometry 



The non-existence of well-defined Cauchy surfaces in regions IV and V of the 
extended space-time means that the canonical theory is impossible to define there. 
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We will therefore not attempt to identify the quantum theory in regions IV and V 
of figure 1. Instead we will ask what are the consequences of spatial difieomorphism 
invariance in the interior. We will see that the wave-functional is required to be 
either vanishing or constant there. 

We will consider the Einstein-Maxwell-dust system described by the action 



S = -±-J d 4 x^[n-F^Fn - ±-J d 4 xy^e(x) [g a/} U a UP + i 

(1.4) 

in the general spherically symmetric space-time 

ds 2 = N 2 dt 2 - L 2 {dr-N r dt) 2 - R 2 dQ 2 , (1.5) 

where N(t,r) and N r (t,r) are respectively the lapse and shift functions, R(t,r) is 
the physical radius or curvature coordinate, e(t,r) is the density of the collapsing 
dust in its proper frame, 1Z is the scalar curvature and U a are the components of 
the dust proper velocity. 

This paper is organized as follows. In section II we summarize the general 
canonical formalism for spherically symmetric space-times, stating the canonical 
form of the action, the appropriate fall-off conditions to be imposed on the canoni- 
cal variables at infinity and the boundary terms. In section III we recast the action 
in terms of a new chart composed of the mass, the curvature coordinate, the dust 
proper time and their conjugate momenta. Sections II and III will closely follow 
Kuchar's original reasoning, 10 which will be adapted to suit the charged black hole 
geometry. In section IV we obtain and solve the Wheeler-DeWitt equation, subject 
to spatial difieomorphism invariance, for the Reissner-Nordstrom black hole, thus 
recovering the spectrum of the black hole. In section V we recast our solution in 
a suitable form and obtain the total entropy and the Hawking temperature in the 
microcanonical ensemble. We conclude in section VI with a few comments on the 
results obtained and the assumptions that went into their making. 

II. Hamiltonian Reduction. 
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The gravitational part of the action in (1.4) was recast by Kuchaf 10 into the 



form 



S 9 



dt l dr 



P L L + P R R - NH9-N r H?\ + S| E , 



(2-1) 



where 



R 

Fi = N 



-R + N r R' 



P " = N 



-LR -LR+ (N r LR)' 
are the momenta conjugate to L and R respectively, and 



(2-2) 



H 9 = 



L 

2 



R 



'2 



2L 



+ 



RR' 



(2.3) 



PlPr LPl 
R 2R 2 

H? = R'P R -LP' L . 

Again, S|| s is a surface term that is required to cancel unwanted boundary terms 
in the variations of the canonical variables and must be determined after specifying 
reasonable fall-off conditions on the canonical variables and the Lagrange multipli- 
ers of the theory. Kuchar's fall-off conditions are well suited to the exterior of the 
maximally extended Reissner-Nordstrom geometry and we shall adopt them here. 
They read 



L(t,r) 


= 1 + M±(t)\r\ 




R(t,r) 


= \r\ + 0°°(\r\- 




P L (t,r) 


= e>°°(r- e ) 




P R (t,r) 


= O^flrr 1 - 6 ) 




N(t,r) 


= N±(t) + C°° 


:\rn 


N r {t,r) 


= 0°°(|rr e ) 





(2-4) 



and, with them, the boundary action required is easily seen to be 

= - J dt[N + (t)M+(t) + N-(t)M-(t)]. 



(2-5) 



Kuchar has emphasized that N±(t) must be considered as prescribed functions of 
the label time coordinate, otherwise a variation of the total action would also lead 
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to the conclusion that the energy of the system at infinity is exactly zero. The 
fact that the N±(t) are prescribed functions will be exploited below to set the 
parametrization clocks at infinity. 

A straightforward canonical reduction of the electromagnetic term in (1.4) with 
the ansatz in (1.5) gives 

S em = J dtjdr [p A A r - NH em - N r H e r m - <f)P A j + Sf%, (2.6) 



where 

R 2 



Pa = 



4 (2.7) 



NL 

is conjugate to A r , 

T P 2 

Tjem _ _|_ J ^ ± A 

" + 2i? 2 (2.8) 
H e r m = -A r P' A , 

and we have defined 

<j>(t,r) = - A t (t,r) + N r (t,r)A r (t,r). (2.9) 

If we adopt the following fall-off conditions 

A r {t,r) = O^dr]- 1 - 6 ) 

P A (t,r) = Q±(t) + C°°(|r|- e ) (2.10) 
0(t,r) = <f>±(t) + 0°°(|r|- e ), 

then the electromagnetic surface term is of the form 



-Jdt - <j>-Q-]. (2.11) 



Q± will turn out to be the electric charge. Once again, to avoid a neutral solution 
(Q = 0) we must treat 4>±{t) as prescribed functions of the label time coordinate. A 
gauge choice that is consistent with the Reissner- Nordstrom solution is 4>±(t) = 
and so this term vanishes. 
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Let us now consider the dust action in (1.4). Dust is described by eight space- 
time scalars, e, r, Z k and W k (k £ {1, 2, 3}). The physical interpretation of these 
variables which follows from an analysis of the equations of motion were given 
in ref. [12] and will be summarized here for completeness, r is the proper time 
measured along particle flow lines, Z k are the comoving coordinates of the dust, 

are the spatial components of the four velocity in the dust frame, and e is 
the dust proper energy density. All these scalars are assumed to be functions of 
the space-time coordinates. In particular, the four variables, Z K = (r,Z k ), are 
independent functions, det|Z^ M | ^ 0, and the four- velocity of the dust particles 
may be defined by its decomposition in the cobasis Z K ^ by 

Ult = -r, + W k Z k ^. (2.12) 

In the spherically symmetric geometry described by (1.5), the dust action may be 
cast into the form 



S d = dt dr 



P T i + P k Z k - NH d - N r Hf 



(2-13) 



where 



LR 2 

Pr = —e(r,t)i-U t + N r U r ] ^ 

P k = -w k p 



T 



are the momenta conjugate to the the dust proper time and the frame variables 
respectively, and 



H d = P T \h+-L 



Hf. = -U r P T . 



L 2 (2.15) 



The expression for H d in (2.15) is obtained upon exploiting the fact that e(t,r) is 
a Lagrange multiplier and therefore 5C/5e = 0. 
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Putting the three components of our system together, we have the Hamiltonian 
form of the total action in (1.4) with the ansatz in (1.5). It reads, 



S = J dt J dr[P L L + P R R + P A A r + P T f + P k Z k -NH -N r H r -<pG] + S 9E , 

(2.16) 

where the boundary action, Sqy, = and is given by the right hand side of (2.5). 
The full super-Hamiltonian and supermomentum constraints are given by 



H 



PlPr LP 2 l 



R 



2R 2 



+ ^7 + Pr\H 



L 
2 



'2 



R_ 
~2L 



0, 



+ 



2R? r V L 2 
R'Pr- LP' l - A r P' A -U r P T w 



RR' 



(2.17) 



and the electromagnetic constraint by, 



G 



P' A » 0. 



(2.18) 



The constraints may be further simplified by requiring that the dust be non- 
rotating and that its motion be described with respect to the frame orthogonal 
foliation. Then we may impose the additional constraints 12 P k = 0. When they 
are applied as restrictions on the state functional ^[r, Z,g,A], they imply that the 
state \& does not depend on the frame variables Z k . The Hilbert space is then 
composed of state functionals, ^[t, g, A], and the quantum theory is described by 
imposing the classical constraints as operator conditions on them, 



H(r,g,A)*[T,g,A] = = H r (r, g, A)*[t, g, A], 



(2.19) 



after a suitable operator ordering has been found. Using (2.12) with W k = 0, the 
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classical constraints take the form 

\P L P R LP 2 
[ R 2R 2 

+ ^ + P T f^ - 0, (2.20) 
H r = R'Pr- LP' l - A r P' A + r'P T « 
G = P' A ~ 0. 

Not only are they not decoupled, making them difficult to solve, but the phase 
space variables are not transparent and "natural" to the black hole problem. The 
spatial hypersurfaces must eventually be embedded in the metric given by (1.1) 
and (1.2) and this line element is completely characterized by the mass, M, and 
the charge, Q, of the black hole. Kuchaf 10 has shown both how these quantities are 
determined by the canonical data as well as how the hypersurface embedding in 
the space-time may be deduced from the values of the phase space coordinates at 
any point. This leads to a reformulation of the constraints (2.20) in terms of more 
transparent variables, the mass, the charge, the physical radius, the dust proper 
time and their conjugate momenta which we describe in the following section. 



L 
2 



'2 



R 
2L 



RR' 



III. New Variables and New Constraints. 

Substituting the foliation, T = T(t,r), R = R(t,r), into the line element given 
in (1.1) and comparing it with the ADM form of the same, i.e., (1.5), we find 

L 2 = -FT 2 + F~ l R 2 
L 2 N r = —FT'T + F~ l R'R (3.1) 
N 2 -L 2 N r2 = FT 2 - F~ l R 2 . 

These identities may be easily solved for the lapse and shift functions, giving 

R'f - T'R 

yZ-FT" 2 + F-^W (32) 
-FT'T + F- X R'R 
-FT' 2 + F~ l R' 2 ' 
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Defined with the positive square- root, the lapse function is positive in all three 
regions of interest with label time going to the future. Substituting the expressions 
(3.2) into (2.2), one finds 

T> = (3.3) 
which can be inserted into the expression for L 2 in (3.1) to give 

d'2 p2 

F(R) = — - ^. (3.4) 

Let us, for the present, work with the mass function M, which we define by F(R) = 
1 — 2M/R. Comparing this with (1.2) it is clear that M must be related to the 
mass and charge of the black hole by 

M = M - ^, (3.5) 

if we are to recover the Reissner-Nordstrom black hole. M as determined by the 
canonical data is 

R 



M 

2 



R^ Pi 
L 2 + R 2 



(3.6) 



Note that M(r) is a local function of the canonical data as is T'(r). If we now 
proceed to compute the Poisson brackets between M(r) and T'(r) from the fun- 
damental Poisson brackets implied by the Liouville form in (2.16) we will see that 
—T'(r) can be interpreted as the momentum conjugate to M(r). Henceforth we 
will refer to it as Pm(j'), i-G., 

TP 

Pm = -T> = (3.7) 

Now, while the pair {M, Pm} has vanishing Poisson brackets with R, it does not 
have vanishing Poisson brackets with Pr. We would like to find a transformation 
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that takes the chart {R, Pr, L, Pl, A r , Pa, t, P t } to a new chart explicitly involv- 
ing the mass and charge of the system, {R, Pr, M, Pm,Q, Pq,t, P r }, subject to 
(3.6), (3.7) and Q(r) = Pa(j'), A r (r) = —Pq{t). The last two relations consti- 
tute an elementary exchange of coordinates and momenta. The four conditions 
provide sufficient information to obtain a generating functional for the canonical 
transformation, which can be given in terms of the original phase space coordinates 
as 



F[R,P R ,L,P L ,A r ,P A ] = Jdr 



A r P A + LP L + — In 



RR' l RR'-LP Li 



, (3.8) 



RR' + LP L 

With the help of (3.6) and (3.7) we compute Pr(t) directly from (3.8); it is 

P R (r) = P R-^-^-^\-( RR ')( LP L)'-(R^)\ LP L)\- (3-9) 

The fall-off conditions can be applied to show that the generating functional in (3.8) 
is well defined near infinity. It can also be shown to stay finite at the horizons. 
On the other hand, the transformation from the old chart to the new is invertible 
everywhere except at the horizons. 

Thus, following Kuchar's reasoning for the Schwarzschild black hole 10 , we have 
introduced the mass and charge as dynamical variables on the phase space. We 
shall now re-express our constraints in (2.20) in terms of the new chart. Again, 
Kuchar has pointed the way: it follows from expression (3.6) for M(r) and expres- 
sions (2.3) that 

M' — — ^-H 9 - ^rHf. (3.10) 
L RL 

This allows us to write the gravitational part of the super-Hamiltonian and super- 
momentum constraints in terms of the new variables. Using the transformations, 
it is easy to see that (2.3) becomes 



H? = R'P R + M'P M 



H 9 = 



M'F~ l R' + FP m Pr 



(3.11) 
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where, expressed in terms of the new variables, 



F = 1 - ™ 

R (3.12) 

L 2 = -FPfa + F^R 2 . 

Furthermore, the constraints H « 0, H r « and G « in (2.20) together imply 
that 

M< M $g + P T */FT^ + ^, (3,3) 



showing that 



M « M - ^, (3.14) 



where M' is given by the last two terms on the right hand side of (3.13). This 
compares with equation (3.5) and shows that to recover the Reissner-Nordstrom 
black hole we must further impose homogeneity via the constraint M'(r) = 0. It is 
not surprising that this condition has to be enforced by hand and does not follow 
directly from the constraints. In introducing non- rotating dust, we have introduced 
an extra degree of freedom in the theory. Thus the problem, as it has been set 
up here, actually describes the more general problem of gravitational collapse of 
inhomogeneous dust before the constraint M'(r) = is imposed. The black hole 
is treated as a special case of the general problem. 

Let us now turn to the boundary term on the right hand side of (2.5). As 
mentioned, N±(t) must be treated as prescribed functions of t. The freedom in 
choosing this function can be combined with the freedom we have of setting the 
dust proper time at infinity to correspond to the parametrization clocks there. The 
lapse function is the rate of change of the proper time with the coordinate time at 
infinity, so we set N±(t) = ±f±(t) to write 



- J dt [M + i+ - M—f. 



Sqy, = - dt [M + i+ - M_f_]. (3.15) 
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It is linear in the time derivatives, f±, and defines a one form, 



[M + 5r+ - M_5r_] = -J dr{M5r)' 



-J dr[M'5r -r'SM + 5(Mr)]. 



(3.16) 



The first two terms on the right hand side may be absorbed into the Liouville 
form of the hypersurface action (they modify the canonical momenta) and the last 
term is an exact form which can be dropped. The action is expressed entirely as a 
hypersurface action. Defining P T = P T — M' and Pm = Pm + r '> we ma Y write it 

as 

S = J dt J dr[P R R + P M M + P Q Q + P T r - NH - N r H r - (3.17) 
where the constraints in the new chart reads 



H 



M'F- X B! + F(P M - r')P R 



+ 



2R 2 



+ {P T + M')\l + 



L 2 



G 



R'P R + M'P M + Q'Pq + r'P T 



Q' 



o, 



(3.18) 



and where L 2 is given in (3.12). A final point transformation, 



R 

Pm 



R, 
Pm, 



M 



Pr 



Q 2 

M + — 
2R 

= Pr + 



Q 2 Pm 
2R 2 ' 



(3.19) 



suggested by (3.14), would express the action above in terms of the true mass, M. 
Furthermore, the Poisson brackets of the new phase space coordinates with the 
constraints would yield the canonical equations of motion in terms of them. 
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If the supermomentum constraint is now used to eliminate Pm m the super- 
Hamiltonian, the latter constraint takes a relatively simple form 



(P T + M'f + F 



2^n 2 



2R 2 



M 



"2 



o, 



(3.20) 



and specializing to the Reissner-Nordstrom black hole by allowing only the homo- 
geneous mode to survive, M'(r) = 0, gives the final form 



2dP 2 



Pr + 



Q l R 
2R 2 



+ F 



2^n 2 



Pr 



Q z r 
2R 2 



(3.21) 



It remains to impose the constraint as an operator equation on the state functional 
ty[r,R,Q] and solve the resulting functional Schroedinger equation. The solutions 
must also obey spatial diffeomorphism invariance and this is imposed by the second 
constraint in (3.18). Acceptable solutions are the topic of the next section. 

IV. Quantization. 

The quantum state of the black hole on a hypersurface r(r), R(r) at a label 
time t is described by a state functional ^[r, R, Q] over the configuration space, 
with the canonical momenta acting upon it as functional differential operators. 
The configuration space can be seen to admit the inverse metric 7 a& = diag(l, F), 
7a?, = diag(l, l/F). The metric is positive definite in the external (and internal) 
regions, and indefinite in the dynamical region. Furthermore, it is a flat metric 
and can be brought to manifestly flat form by a coordinate transformation. In the 
external region, this transformation takes the form 



R* — 



dR 



RdR 



y/F{K) J ^R 2 - 2MR + Q 2 
= R\jF(R) + M\n\R-M + Ry/F(R)\ 



(4.1) 
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and in the dynamical region it is 

dR 



R* 



/ 



/ 



RdR 



V^FjR) 
-Ry/-F(R) -Msin- 1 



^2MR -R 2 -Q 2 
M-R 



a/m 2 - q 2 



(4.2) 



These definitions may be modified by additive constants so as to turn i?* into a 
continuous variable. The classical constraints can now be expressed in terms of the 
momentum conjugate to i?* in each region. They are 



i^P* + M'P M + Q'Pq + r'P 7 



Pr + 



Q 2 R' 
2R 2 



± 



Q 2 ^±F\R)r' 
2R 2 





0, 



(4.3) 



where the ± signs refer to the exterior and the dynamical region respectively, and 
must be imposed as operator constraints on the wave functional. In the new, 
manifestly flat, configuration space, we represent the momenta by the functional 
derivatives 

f T (r) = -i 
%(r) = -i 
and consider the solutions of the equation 



5r(r) 
5 

SRJr) 



(4.4) 



5 T + i 



Q Z R 
2R 2 



± 



5* —i 



Q 2 ^±FJR)t> 
2R 2 



V[t,R*,Q} = 



(4.5) 



that respect the diffeomorphism constraint, 



{R% + r'5 T }^[r,R„Q] = 0, 



(4.6) 



where we have already imposed the two requirements, Q'^ = and M'^ = 0, the 
first of which is simply the electromagnetic constraint and the need for the second 
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arises, as we have mentioned earlier, from the fact that we have added a degree of 
freedom to the system in the form of dust, which must be constrained in order to 
describe the Reissner-Nordstrom black hole. Equations (4.5) and (4.6) define the 
quantum theory whose Hilbert space is H := £ 2 (R, dR*) with inner products 

oo 

(*i,* 2 ) = J dRJ&\v 2 (4.7) 

M In a/M 2 -Q 2 

in the exterior region and 



2 

(*i,* 2 ) = J dR^ 2 (4. 



7fM 

2 



in the dynamical region. 

Now it is readily verified, by taking functional derivatives, that any solution of 
equation (4.5) can be written in the form 



*o[t,R*\, (4-9) 



^[t,R*,Q] = exp 

where \l/ is a solution of the free equation 

(6 2 T ±6t)* [T,R*} = 0. (4.10) 

We have, once again, used the ± signs to refer to the exterior and the dynamical 
regions respectively. Not surprisingly, (4.10) was obtained for the Schwarzschild 
black hole in refs. [8,9]. Two features of this equation are worthy of mention. For 
one, a consequence of the signature change in the configuration space metric as 
we move from the external region to the dynamical region is that the "equation of 
motion" goes from being elliptic to hyperbolic. For another, (4.10) is decoupled 
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implying that it may be solved independently at each point, labeled by r, of the 
spatial hypersurface. In the exterior, then, the solutions will be exponentially 
decaying, but they will oscillate in the dynamical region and in each region they 
will have the form 

*o[t,R*\ = H^ [r(r),R4r)}. (4.11) 

r 

Thus we may write a general solution of (4.5) as 

H^ [r(r),R4r)}. (4.12) 

r 

Again, the solution will obey spatial diffeomorphism invariance if and only if 
V [t,R*\ obeys (4.6). 

Consider the solutions in the exterior (F > 0). The general positive energy 
solution that is well behaved in the entire range of R* is given by 

^ [r(r),R,(r)] = c(M, Q) exp [-iE(r(r) - iR*(r))] , (4.13) 

where c(M, Q) is a mass and charge dependent constant. However, upon consid- 
ering the action of the spatial diffeomorphism invariance constraint on ^ofr, R*] 
defined in (4.11), we conclude that Eir' — iR'^)ty = 0. For positive energy solu- 
tions, this condition is met by r' = = R'^, by E = or by c(M, Q) = 0. However, 
R* = implies from (4.1) that R' — and this is unacceptable in the exterior 
because L 2 given in (3.12) is required to be positive definite. We could take E = 
and ^ would be a constant, c(M, Q), but this solution would not be C? according 
to (4.7). We conclude that c(M,Q) must vanish and consequently that the wave 
functional, ^ [t,R^], is identically zero in the exterior of the black hole. 

In the dynamical region, F < 0, the situation is completely different. The 
"equation of motion" is hyperbolic and the solutions are oscillatory. The general 



^[t,R*,Q] = exp 



i J dr ^ 



Q 2 

2R(r) 



r(r) 
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positive energy solution is now given by 



^ [r(r),R^r)] = a(M, Q) e -iE[r{r)+R,{r)] + ft ( M> Q^E[r(r)-R t (r)] > ( 414 ) 

where a(M, Q) and 6(M, Q) are mass and charge dependent constants. Because 
the wave functional in the exterior vanishes, ^ f [r(r), R*(r)] in (4.14) must vanish 
on the outer horizon, R* = nM/2. This gives 

b(M,Q) = -a(M,Q)e- inEM (4.15) 



and 



V [r(r),R*(r)] = a(M,Q) \ e ~ iE ^ +R *^ - e ~^EM e -iE[r(r)-R,(r)}] (416) 

Further boundary conditions come from matching the wave-functional, ^[r, R*] 
above, with its counterpart in the interior on the inner horizon. As we have men- 
tioned in the introduction, the inner horizon is a Cauchy horizon for spatial sec- 
tions. It is impossible as far as we know to define the dynamics consistently here 
because the Cauchy data must be supplemented by boundary conditions on the 
singularities which are impossible to give. We would expect that, since there is no 
dynamics within the inner horizon and since the diffeomorphism constraint must 
hold there, the wave function will be constant or equal to zero everywhere. If we 
take it to be identically zero and require that ^f [r, i?*] vanish on the inner hori- 
zon as well, one sees that the energy is quantized according to the simple relation 
EM = n, where n is a non-negative integer. 

Yet, a vanishing wave- functional in the interior may not be the only possibil- 
ity. Even though we can make no statements about the evolution equations, it is 
reasonable to expect that spatial diffeomorphism invariance of the wave-functional 
in the interior is respected according to (4.6). In the privileged foliation defined by 
the dust proper time, r' — but R' ^ 0. In fact, R' can never be zero in the static 
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regions. It follows from (4.6) that the wave- functional on constant proper time 
hypersurfaces must be constant, i.e., it may depend only on r. One easily verifies, 
for example, that the charge dependent phase of (4.9) vanishes on every leaf of 
this foliation. In order for the wave-functional in the interior to be consistently 
matched with the wave-functional in the dynamical region it must have the form 

¥ nt [r,R*} = l[d(M,Q)e- tET . (4.17) 

r 

This, of course, does not respect the super-Hamiltonian constraint in (4.5), but 
we do not wish to impose it in the interior. Matching the wave-functionals in the 
interior and in the dynamical region one finds that the energy is quantized in half 
integer units of the Planck mass (as opposed to integer units when d(M, Q) = 0) 
according to EM = (n + 1/2). The ground state energy is non-zero excluding, as 
we will see, the extremal solution as well as the Minkowski vacuum. Moreover, it 
does not reduce to the spectrum of the Schwarzschild black hole 8 ' 9 in the Q — > 
limit. It will not change our counting of the black hole states and we will not 
pursue this case further in this paper. 

We must now address the physical meaning of the momentum conjugate to the 
dust proper time in region III, i.e., what is the functional dependence of the proper 
"energy", E, on the mass and charge of the hole? By our coupling to dust, the 
time coordinate is always chosen to coincide with the proper time, r, of a freely 
falling observer. It can be expressed in terms of the curvature coordinates as 

r(T,R) = T + J dR ^ l n ^ R) (4.18) 
in the (static) regions I, II, IV and V, and 

T (R) = R* = I dR (4.19) 

in the dynamical region, III. Transformations (4.18) and (4.19) generalize the cor- 
responding transformation for the Schwarzschild black hole given in ref. [14]. A 
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proper time interval in the dynamical region is seen to agree with the proper time 
interval of the asymptotic observer and with his Minkowski coordinate time. Its 
conjugate momentum, or the proper energy, will therefore be the energy contained 
in this region, between the inner and outer horizons. Now, the exterior and interior 
of the space-time admit a time-like Killing vector, so it is natural to think in 
terms of Komar's definition 13 of the energy 

E = — 

S 

where S is a spatial two sphere. For an exterior which is vacuum, for example 
the Schwarzschild black hole, the integral is independent of the surface, S, and 
the definition requires only that ^ is a time like Killing vector. If the surface is 
chosen to lie at infinity, Komar's definition may be applied to all asymptotically 
flat space-times to obtain the total energy. In regions I, II, IV and V, for a two 
sphere S with curvature radius R, equation (4.20) may be thought of as the total 
energy, E, of the system interior to S. The energy contained between the horizons 
would be the difference between the energies interior to the two horizons, i.e., its 
energy interior to the outer horizon minus its energy interior to the inner horizon. 
Therefore, we define the energy E in the dynamical region as 

E = -ll tdi + il*^ (421) 

S+ S- 

where S- and S+ refer to the inner and outer horizons respectively with radii 
R T = Mt \/M 2 - Q 2 . By its definition, E is a conserved quantity and a true 
invariant. It gives E = \/ M 2 — Q 2 = (R + — R-)/2 and reduces appropriately to 
E = M in the limit as Q -> 0. 
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Fig. 3: The extremal Reissner-Nordstrom geometry 

Furthermore, region III disappears for the extremal black hole (see figure 3), as 
does the proper energy associated with it. 

Applying the definition of the energy, E, appropriate to this region we find the 
quantization rule 

My/ M 2 - Q 2 = nM%, (4.22) 

where we have introduced the Planck mass, M p . We have also thus recovered the 
Bekenstein mass spectrum for the Schwarzschild black hole (Q — > 0), which was 
examined separately in refs. [8,9]. Moreover, if A± refer to the areas of the outer 
and inner horizons respectively, the quantization condition reads 

A+ - A- = 167rn^, (4.23) 

i.e., it is the difference between the outer and inner horizon areas that is quantized 
in integer units. 

We must of course ensure that the wave-functionals in this region obey the 
supermomentum constraint (4.6). This is easily verified. Applying (4.6) on the 
wave-functional shows that r'(r) = = R'(r), assuming that E ^ 7^ a(M,Q). 
Both these conditions can be easily met in the dynamical region. Here there is no 
contradiction with the positivity of L 2 because F < and we see that the phase 
in (4.12) disappears leaving only the direct product state in (4.11). The solutions 
can be classified as even parity and odd parity states, just as in the case of the 
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Schwarzschild black hole, 



*(+)[r,i2*] = Yl-^=e- iET ^coa[EIU(r)] EM = (2n + 1) , 

r \ M (4.24) 
*(-)[r,i2*l = lT^=e" i£;r(T ' ) sm[ J E; j R4r)] EM = 2n , 



where we take n 6 N U {0} in keeping with the positive energy requirement and 



E = \/M 2 - Q 2 . Then if we think of the extremal Reissner-Nordstrom black hole 
as the limiting case of the non-extremal black hole, we see that it corresponds to 
n = with a vanishing wave functional. As a consequence, the entropy of the 
extremal black hole vanishes identically. It would appear that the extremal black 
hole either has no dynamics at all or that it cannot be understood as the limit 
Q — > M of the non-extremal black hole and must be treated independently. It 
would be of interest, for example, keeping the exterior fixed, to consider different 
interior solutions that leave the horizon non-singular and contain modes which 
propagate entirely within the horizon. These modes will contribute to the entropy 
of an extremal black hole, giving a non-zero result. This is the approach of ref. [15], 
but how this program can be implemented for the non-extremal case and within 
the context of the canonical theory is not clear to us at present. 



V. The Entropy. 

In the previous section we have seen that the states of the black hole reside 
only in the dynamical region, R_ < R < R + , between the inner and outer horizons 
of the black hole and that the Wheeler-DeWitt equation is decoupled so that the 
wave-functional is expressible as a direct product state 

r 

where each component of the direct product over the label coordinate is either an 
even parity state or an odd parity state in (4.23). We may imagine that a lattice 
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is placed on the spatial hypersurface so that the classically continuous variable r is 
a discrete label. The Wheeler-DeWitt wave functional then represents a collection 
of, say, N (assumed to be finite) decoupled oscillators each determined by the same 
Schroedinger equation and obeying the same boundary conditions. The black hole 
entropy is a consequence of the fact that a knowledge of its total mass and charge is 
not equivalent to a knowledge of the number of ways in which this mass and charge 
is distributed between the components. Each particular distribution corresponds to 
a definite microstate of the black hole. To compute the entropy we must enumerate 
these microstates. 

It is therefore convenient to reformulate the problem by recognizing that the 
wave equation in the dynamical region at each label coordinate is derivable from 
the action 

TTM 

2 

Sr = ~\ J ^X^Y'da^ld^r, (5.2) 

7rM 

2 

where X e (r, /?*), and that the total action has the form 

N N 2 

S = J^Sr = J ^X^^da^d^r. (5.3) 

r =l r=1 _zM 

2 

The boundary conditions are that each \l/ r vanishes at the outer and inner horizons. 
Performing a mode expansion of ^ r and combining both parities, we express the 
contribution of any one of the lattice sites to the total energy of the system in 
terms of pairs of creation and annihilation operators, (a n ,a n ) and ((3 n ,[3 n ), 

[a n , a m ] = n5 nm 

(5.4) 

as follows 

M 2 

= -JlY.( a ^ a n r + PjPn r ). (5.5) 
n r 
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The total energy is the sum over contributions from each of the lattice sites, i.e., 

M 2 N 

r r—1 n r 

and this must be the energy of the dynamical region. We consequently obtain a 
dispersion relation of the form 

M 2 N 

E = wJ2J2^ N nr + w 

r=l n r ,l r ^p, 
M 2 0- M l r 

M ^ M 1 1 

r=l 

where n r , l r and iV nr , A^/ r refer respectively to the level number and the occupation 
number at level number n r (l r ) corresponding to the oscillator at lattice site r. 
Equation (5.7) is yet another way of seeing that it is the product, EM, that is 
quantized in integer units. 

A crude estimate of the number of states may be given as follows. Let p(u r ) be 
the density of levels describing each site r. It will naturally account for the angular 
variables which have been integrated over in the canonical theory and have so far 
played no role. Assuming that the lattice sites are distinguishable, for a generic 
level density and large v r the density of states can be written as 9 



N °? / N \ 

Q(E, M, N) = Y[ / du r p(u r )5 EM - ^ v a I , (5.8) 

r=d V s=l I 



where v Q represents a lower limit on the energy attributable to each oscillator, 
which we take to be zero. The 5— function restricts the limits of the v r integrals in 
the product. Following Carlitz 16 and Frautschi 17 we estimate the r th integral by 

J dv r p{y r \ (5.9) 
o 
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provided that the upper limits, A r , are subject to the condition 

N 

= EM. (5.10) 

r=l 

Then the maximum contribution to Q(E, M, N) is obtained when all of the A r are 
of the order EM/N. This provides an estimate for the integrals in (5.9). Quite 
generally 

Q(E,M,N) = a N f N (0, (5.11) 

where £ = bEM/N and a and b are constants which can be determined from the 
density of states. The precise value of these constants is irrelevant because, as we 
shall see, the functional dependence of the entropy on the mass and charge of the 
black hole is a consequence only of the dispersion relation, EM = Ylr=i v r-> m 
(5.7). So far iV has been introduced by hand as the number of lattice points on a 
typical hypersurface. We determine its value by maximizing the number of states 
with respect to it. Notice that 

— \nSl(E,M,N) = ma + ln/(0 - = (5.12) 

is a homogeneous equation in £ whose solution will yield a value, say £ = a -1 . 
Assuming that it is real, this value determines the number of lattice points as 

N max = abEM (5.13) 

and thus the entropy as 



S(M,Q) = \nQ(E,M,N max ) = 1 M V / M 2 - Q 2 , (5.14) 
where 7 is determined in terms of a, a and b, by inserting (5.13) in (5.11), 

7 = abln [a/(a -1 )] . (5.15) 
The entropy is the difference between the outer horizon area and the inner horizon 
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area, 

and is quantized in integer units. Our result does not coincide with Bekenstein's 
proposal, S A+, except if Q — > 0. 

VI. Discussion. 

We have generalized our earlier study of the energy spectrum and the statistical 
entropy of the Schwarzschild black hole to the charged, Reissner-Nordstrom black 
hole. The canonical reduction differs only slightly from Kuchar's treatment of the 
Schwarzschild black hole and the quantization program is seen to lead to precisely 
the same wave- functional as was obtained in refs. [8,9] for the Schwarzschild black 
hole, when it is taken to vanish in the interior, static region. Our fundamental result 
is that if E is the energy associated with the dynamical region, the quantization 
condition and entropy are given by 

EM = MM 2 D , N eNU{0}, 

(6.1) 

S = jEM. 

where 7 is a constant that can be determined. The key distinction between the 
charged and uncharged black hole is in what we define to be the "energy" associated 
with the dynamical region. In the absence of a canonical choice and motivated by 
physical considerations, we have made what we consider to be a reasonable propo- 
sition: take it to be the energy contained between the inner and the outer horizons 
of the Reissner-Nordstrom solution and calculated as the difference between the 
Komar integrals evaluated on the horizons in the exterior and in the interior, where 
time-like Killing vectors exist. This then leads to the conclusion that it is the dif- 
ference in the horizon areas that describes the entropy and that is quantized in 
integer units. The result does not agree with Bekenstein's hypothesis, but it does 
confirm the existence of an "area quantization" law. 
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It follows that the temperature of the charged black hole, as measured by an 
observer at infinity, does not agree with its semi-classical Hawking temperature. 
Indeed it is lower than the Hawking temperature due to a contribution from the 
inner horizon. One finds, 



_! 7 R 2 + + R 2 _ 

' : \m) Q ~ 2 V /M2-Tg2" 



The first term on the right hand side is the inverse Hawking temperature and is 
attributable to the outer horizon. Likewise, the second term is an inverse "tem- 
perature" that is attributable to the inner horizon. This inner "temperature" is 
negative. If we call these temperatures, respectively, T + and TL, the temperature 
of the charged black hole can be written as 

111 . . 

T = %~ T~J (6 - 3) 

where Tg 1 = (dS±/dM) Q and S± = 7^ ± /16tt. 

The canonical quantization described in this article leads to the amusing picture 
of a black hole as a microcanonical ensemble of oscillators, similar to hadrons whose 
statistical mechanics was studied many years ago by Carlitz 16 and Frautschi 17 . The 
oscillators are imagined to be located rigidly on a lattice with N sites placed on a 
spatial hypersurface, and the optimal number of sites is determined by maximizing 
the density of states. It turns out that this number is proportional to the difference 
in horizon areas, which leads to the expression we derived for the entropy. N is the 
number of dynamical degrees of freedom of the system and this statistical result 
encourages the following interpretation. Each horizon area statistically yields the 
total number of degrees of freedom within it. An outside observer sees the total 
number of degrees of system within R + minus the total number of degrees of 
freedom within R-, which difference is the number of degrees of freedom added 
dynamically. 
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What physical processes are responsible for the "lost" degrees of freedom? 
The temperature and the entropy of the black hole are eventually measured by 
an external, distant observer by placing a "thermometer" in the form of a test 
quantum field in the black hole background. Only those modes of the quantum 
field that do not enter the outer horizon are accessible to the asymptotic observer. 
A transition from a higher energy state to a lower energy state by the emission 
of a single quantum of radiation may occur on any oscillator, including one whose 
degree of freedom is attributable to the inner horizon. However, the quantum field 
in the exterior would not be sensitive to an emission from such an "inner" oscillator 
so such a transition would elicit no response from the "thermometer" . As far as 
the asymptotic observer is concerned, therefore, these degrees of freedom must be 
subtracted from the total. The semi-classical theory does not account for this loss, 
hence the expression (5.16) for the entropy is lower than the semi-classical analysis 
would predict and the temperature associated with the black hole in the canonical 
theory lower than the Hawking temperature. 

An equivalent way of viewing this is to observe that the semi-classical analysis 
does not involve the back reaction to any radiation emitted exterior to the black 
hole. Electrically neutral radiation, for example, would diminish the mass of the 
black hole but leave its charge unchanged. Hence R+ would decrease but R- 
would increase. Semi-classically, the asymptotic observer would be unaware of this 
increase. He would also be unaware of a similar effect induced by the possible 
emission of radiation from the interior of the black hole, which radiation would 
either be absorbed by the singularity or pass on to another branch of the universe. 
Hence the measure of entropy and temperature given by the asymptotic observer 
appears to need adjustment or interpretation. 

While it may be argued, as we have, that the semi-classical analysis is not 
sensitive to the conditions in the interior of the hole and is therefore not reliable 
when multiple horizons are involved, it is also possible that the semi-classical result 
is correct and that the boundary conditions we have imposed are too restrictive. 
We have treated the interior as the analytical extension of the exterior, but it has 
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been suggested , for example, within the string theory framework and for extreme 
black holes, that the internal region must be replaced by solutions which contain 
modes propagating entirely within it and matched smoothly to the exterior. At 
the time of writing, we do not know how such a program would be implemented 
for non-extremal black holes within the context of the canonical theory. 
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